arXiv:1507.04554vl [physics.soc-ph] 16Jul2015 


Traffic-driven epidemic spreading in correlated networks 

Han-Xin Yang, 1 El Ming Tang, 2 0 and Ying-Cheng Lai 3 

1 Department of Physics, Fuzhou University, Fuzhou 350108, China 
2 Web Sciences Center, University of Electronic Science and Technology of China, Chengdu 610051, China 
3 School of Electrical, Computer and Energy Engineering, Arizona State University, AZ 85287, USA 

(Dated: July 19, 2015) 

In spite of the extensive previous efforts on traffic dynamics and epidemic spreading in complex networks, 
the problem of traffic-driven epidemic spreading on correlated networks has not been addressed. Interestingly, 
we find that the epidemic threshold, a fundamental quantity underlying the spreading dynamics, exhibits a 
non-monotonic behavior in that it can be minimized for some critical value of the assortativity coefficient, a 
parameter characterizing the network correlation. To understand this phenomenon, we use the degree-based 
mean-field theory to calculate the traffic-driven epidemic threshold for correlated networks. The theory pre¬ 
dicts that the threshold is inversely proportional to the packet-generation rate and the largest eigenvalue of the 
betweenness matrix. We obtain consistency between theory and numerics. Our results may provide insights 
into the important problem of controlling/harnessing real-world epidemic spreading dynamics driven by traffic 
flows. 

PACS numbers: 89.75.Hc, 05.70.Ln, 05.60.-k 


I. INTRODUCTION 

Both epidemic spreading CHU and traffic transporta¬ 
tion 11 2 h 1811 are two types of fundamental dynamical pro¬ 
cesses on complex networks. The past fifteen years have wit¬ 
nessed a great deal of effort in exploring and understanding 
these processes. Intuitively, spreading should be greatly af¬ 
fected by the transportation dynamics as the traffic flows de¬ 
termine the paths along which a virus or information can prop¬ 
agate. Nevertheless, for a long time, the two types of dynami¬ 
cal processes were studied separately. The first series of works 
to incorporate traffic dynamics into epidemic spreading were 
done using the metapopulation model lfl9l - {25tl . In particular, 
consider a spatially extended system of a large number of in¬ 
dividuals. In a metapopulation model, all the individuals are 
grouped into a number of spatially structured subpopulations 
interacting with each other in a manner that can be described 
as a network. The links in the network represent the traveling 
paths of individuals across different subpopulations. An in¬ 
fected individual can infect other individuals in the same sub¬ 
population. The metapopulation model has been used to sim¬ 
ulate the spread of human and animal diseases (such as S ARS 
and H1N1) among different cities. In 2009, Meloni et al. pro¬ 
posed another traffic-driven epidemic spreading model [26], 
where each node corresponds to a router (e.g., in a computer 
network) and the virus can spread among nodes through the 
transmission of packets. A susceptible node will be infected 
with certain probability every time it receives a packet from 
an infected neighboring node. This model is quite suitable for 
studying computer virus propagation. 

The model of Meloni et al. [26] has become a prototypi¬ 
cal framework to address a variety of issues in traffic-driven 
epidemic spreading. For traffic flow dynamics on complex 
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networks, an essential ingredient is the routing strategy. In 
the model, one can demonstrate that epidemic spreading can 
be modulated or controlled through a proper choice of the lo¬ 
cal routing strategy [[270 or global routing protocol l28t] . For 
a fixed routing strategy, the network structure can also affect 
the spreading dynamics. For example, by increasing the aver¬ 
age network connectivity, one can effectively delay epidemic 
outbreaks Furthermore, the epidemic threshold can be 
enhanced by deliberately removing a subset of links associ¬ 
ated with the large-degree nodes or the edges with the largest 
algorithmic betweenness [30], 

In previous studies of traffic-driven epidemic spreading, the 
networks were assumed to be uncorrelated. In such a network, 
the average degree of the neighboring nodes is uniform across 
nodes. However, many real-world networks display various 
degrees of mixing patterns | f3lil . Qualitatively, a network is as- 
sortatively (disassortatively) mixing if high-degree nodes tend 
to connect with high-degree (low-degree) nodes. Quantita¬ 
tively, degree mixing can be characterized by the assortativity 
coefficient & 

= M ~ l E i Jih - [M- 1 \{ji + k)F 

C M-1 + kf) - [M-1 £i + h)} 2 ’ 

where j, : , ki are the degrees of the nodes at the ends of the 
vth edge, M is the number of edges in the network, and i = 
1,..., M. For standard network models such as the Erdos- 
Renyi random graphs {32,] and the Barabasi-Albert scale-free 
networks [33], the assortativity coefficient c is zero, indicat¬ 
ing complete lack of degree correlation. For many social net¬ 
works, the assortativity coefficient c is positive. However, 
technological and biological networks tend to be disassorta- 
tive with negative values of c. 

In this paper, we study quantitatively how degree mixing, 
the most pronounced feature of correlated networks, affects 
traffic-driven epidemic spreading. Interestingly, we find a 
non-monotonic behavior in that the threshold can be mini¬ 
mized for some critical value of the assortativity coefficient. 
To understand this phenomenon, we employ a degree-based 
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mean-field theory, which allows us to calculate the traffic- 
driven epidemic threshold in correlated networks. We find 
that the threshold is inversely proportional to the packet- 
generation rate - a key parameter characterizing the traffic 
dynamics, as well as the largest eigenvalue of the between¬ 
ness matrix. Moreover, the epidemic threshold is negatively 
correlated with the maximum algorithmic betweenness. Our 
theory predicts the existence of a critical value of the assorta- 
tivity coefficient at which the epidemic threshold reaches its 
minimum, which agrees with the results from direct numerical 
simulations. 


H. CORRELATED SCALE-FREE NETWORK AND 

TRAFFIC-DRIVEN EPIDEMIC SPREADING MODELS 

Construction of scale-free networks with tunable assorta- 
tivity coefficient. We first generate an ensemble of scale- 
free networks according to the uncorrelated configuration 
model lf34ll . The steps are: (i) assign to each node i, in a set of 
N initially disconnected nodes, a number ki of stubs, where 
ki is drawn from the probability distribution P(k) ~ A; -7 and 
subject to the constraints k m i n < ki < s/N (here we set 
kmin = 2 and N = 4900) and JT fc, being even, and (ii) 
generate a network by randomly choosing stubs and connect¬ 
ing them to form edges, taking into account the preassigned 
degrees and avoiding multiple and self-connections. We then 
use the algorithm in Ref. [35] to obtain networks with desired 
degree mixing patterns. Specifically, in order to obtain an as- 
sortative network, at each step we randomly choose two dif¬ 
ferent edges with four different ends, and then purposefully 
swap the two edges by linking the vertices with higher de¬ 
grees and lower degrees, respectively. Repeating this proce¬ 
dure while forbidding multiple connections and disconnected 
components, we generate a network with certain degree as- 
sortativity without altering the degree distribution of the orig¬ 
inal network. Through the opposite operation that one edge 
links the highest and the lowest nodes and the other edge con¬ 
nects the two remaining nodes, we can obtain networks with 
disassortative mixing. When the assortativity coefficient of 
a network reaches some expected value, we stop the edge¬ 
swapping process. 

Traffic-driven epidemic spreading model. Following the 
pioneering work of Meloni et al. a, we incorporate 
the traffic dynamics into the standard susceptible-infected- 
susceptible epidemic spreading model lf36ll . as follows. At 
each time step, R new packets are generated with randomly 
chosen source and destination nodes, and each node i can de¬ 
liver at most Di packets toward their destinations. (For sim¬ 
plicity, we set Di to be infinite.) Packets are forwarded ac¬ 
cording to a given routing algorithm (we use the shortest-path 
routing scheme). The queue length of each node is assumed to 
be unlimited, to which the first-in-first-out principle is applied. 
Each newly generated packet is placed at the end of the queue 
of its source node. Once a packet reaches its destination, it 
is removed from the system. A node can be in two discrete 
states: susceptible or infected. After a transient time, the to¬ 
tal number of delivered packets at each time step will reach a 



FIG. 1: (Color online) For scale-free networks of degree exponent 
7 = 3 with different values of c, the average degree of the nearest 
neighbors, (k nn (k)), as a function of the degree k. Each data point 
is the result of averaging over 100 independent network realizations. 
For assortative networks (c > 0), ( k n n(k )) is an increasing func¬ 
tion of k (for small k values). For disassortative networks (c < 0), 
( k n n{k )) is a decreasing function of k. For uncorrelated networks 
(c = 0 ), ( knn{k )) is independent of k. 


steady value. Subsequently, an initial fraction of nodes p 0 is 
set to be infected (we choose po = 0.1 in numerical experi¬ 
ments). Differing from the conventional spreading dynamics 
where there is a certain probability that a node can be infected 
if one or more of its neighboring nodes are infected, here the 
infection spreads in the network through packet exchanges. In 
particular, all packets queuing at an infected node are infected, 
while all packets at a susceptible node are uninfected. A sus¬ 
ceptible node has the probability /3 of being infected every 
time it receives an infected packet from any infected neigh¬ 
boring nodes. With probability 1 — /3, the virus in an infected 
packet will be cleaned by some anti-virus software (say, in 
a computer network) in the susceptible node. The infected 
nodes are recovered at the rate p (we set //, = 1). 


III. NUMERICAL RESULTS AND MEAN-FIELD THEORY 

Numerical results of the dependence of the epidemic thresh¬ 
old on the assortativity coefficient. A useful quantity to char¬ 
acterize correlated networks is the average degree of the near¬ 
est neighbors, denoted as ( k nn {k )). Figure Q] shows ( k nn {k )) 
as a function of degree k for different values of the assorta¬ 
tivity coefficient c. From Fig. Q] one can see that (k nn (k)) 
is almost independent of k when c = 0. For assortative (dis¬ 
assortative) networks, ( k nn (k )) is an increasing (decreasing) 
function of k, indicating that large-degree nodes tend to con¬ 
nect with large(small)-degree nodes. 

A fundamental quantity in SIS dynamics is the epidemic 
threshold ffi, below which the epidemic dies off asymptot¬ 
ically [37, [38]. Figure U shows /3 C versus the assortativity 
coefficient c for different values of the scale-free degree ex- 
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ponent 7 . We see a non-monotonic behavior. Especially, for 
each value of 7 simulated, there exists a critical value of c, de¬ 
noted as c wor , at which the epidemic threshold is minimized. 
The inset of Fig0shows that c wor decreases from 0.3 to 0.15 
as 7 increases from 2.1 to 4. 

Mean-field analysis of epidemic threshold for correlated 
networks. We seek to understand the counterintuitive behav¬ 
ior exemplified in Fig. [2] The base of our analysis is the stan¬ 
dard degree-based mean-field theory [39] for complex net¬ 
works. First, we write the dynamical rate equations for the 
traffic-driven SIS model as 

dpk ^ = -Pk(t)+(3kR[l-p k (t)]^2b k '^ k P(k \k)p k -{t), 

k' 

( 2 ) 

where pk is the density of the infected nodes of degree k, 
P(k | k) is the conditional probability that a node of degree 
k is connected to a node of degree k , and b k '^ k is the algo¬ 
rithmic betweenness (40] of a directed edge from a node of 
degree A; to a node of degree k. The first term in Eq. (0 is the 
recovery rate of the infected nodes, and the second term takes 
into account the probability that a node of degree k belongs 
to the susceptible class of fraction 1 — pk(t) and gets infected 
via packets from some infected nodes. The algorithmic be¬ 
tweenness of a directed edge, bi^j, is the number of packets 
passing from node i to its neighboring node j each time if the 
packet-generation rate is R = 1. For an undirected network, 
we have bi^j = For the shortest-path routing protocol, 
the algorithmic betweenness of a directed edge bi^j can be 
calculated through 


bi ^ 3 N(N - 1) ' 

V ’ m^n 


3) 


(3) 


where a m ->-n is the total number of the shortest paths going 
from node m to n, and cr m ^ n (i —> j) is the number of short¬ 
est paths going from m to n and passing from node i to its 
neighboring node j. The algorithmic betweenness of a node i 
is given by b, = where the summation runs over 

the nearest neighbor set f \ of node i. 

For small values of pk, we can linearize Eq. 0 to obtain 

dpk(t) _ ^ T ti\ iA\ 

^ — 2^ ^ kk ' P k ' O’ (4) 

k' 

where the Jacobian matrix element is 

L kk' = -<W + PkRP{k\k)b k '^ k , (5) 

and S kk > is the Kronecker delta function. The solution pk = 0 
will be unstable if there exists at least one positive eigen¬ 
value of the Jacobian matrix L. The endemic phase (pk > 0) 
will thus take place for — 1 + /3 RA m > 0, where A m is the 
largest eigenvalue of the betweenness matrix C with its ele¬ 
ments given by 

C kk‘ = kP ( k '\ k ) b k'^k- ( 6 ) 

The element C kk > gives the number of packets that a node of 
degree k receives from all neighboring nodes of degree k at 



FIG. 2: (Color online) For scale-free networks of different values 
of degree exponent 7 , epidemic threshold as a function of the 
assortativity coefficient c. The solid, dashed, and dotted curves are 
theoretical predictions for the cases 7 = 2.5, 3, and 3.5, respectively. 
The inset shows the critical value c WO r as a function of 7 . The packet- 
generation rate is R = 1500. Each data point is result of averaging 
over 100 random network realizations. 
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FIG. 3: (Color online) For scale-free networks of degree exponent 
7 = 3 and different values of the assortativity coefficient c, (a) al¬ 
gorithmic betweenness bk and (b) density of infected nodes pk as a 
function of the degree k. For all cases, the packet-generation rate is 
R = 1500 and the density of the infected nodes in the whole popu¬ 
lation is p ~ 0.04. Each data point results from averaging over 100 
random network realizations. 
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FIG. 4: For scale-free networks of degree exponent 7 = 3 , dependence of the infection probability pi on the algorithmic betweenness bi 
for different nodes. The values of the assortativity coefficients c are -0.2, 0 and 0.6 for the left, middle and right panels respectively. Other 
parameters are: packet-generation rate R = 1500 and density of infected nodes p « 0.04. 


each time step when the packet-generation rate R = 1. Since 
C is non-negative and irreducible, the Perron-Frobenius theo¬ 
rem El stipulates that its largest eigenvalue is real and posi¬ 
tive. The endemic state then occurs for 


/3 > Pc 


1 

f?A m ' 


(7) 


For uncorrelated networks, we can use two approximations: 
( i ) P(k |fc) oc k P(k ), and ( ii ) b k >^ k oc bkb k > /kk , where 
the latter allows us to identify the key parameters of the un¬ 
derlying dynamics. Under these approximations, the element 
of the betweenness matrix can be reduced to 


CH, « b k .b k P{k')/{b), 

which has a unique non-zero eigenvector with the eigenvalue 
A“ c = (b 2 ) / (b). The epidemic threshold for uncorrelated net¬ 
works can then be obtained as 


]UC _ (^) 

c “ R(b 2 ) ’ 


( 8 ) 


where b represents the algorithmic betweenness of a node. 
This result was obtained previously by a slightly different ver¬ 
sion of the mean-held theory (heterogeneous mean-held the¬ 
ory) 0 . 

To validate our theoretical results [Eq. 0 ], we compare 
them with numerically calculated dependence of /3 C on the as¬ 
sortativity coefficients c. As shown in Fig. [2j the theoretical 
and numerical results are consistent. 

How degree mixing affects the algorithmic betweenness 
and the infection probability for nodes of different degrees? 
The nodal algorithmic betweenness within the degree class k 
is dehned as bk = T,ibi/N k , where N k is the number of nodes 
within each degree class k and the summation runs over all 


nodes of degree k. Figure Oa) shows the dependence of b k 
on degree k. We see that for dissortative (e.g., c = —0.2), 
uncorrelated (c = 0), and weakly assortative (e.g., c = 0.25) 
networks, bk increases with k , and the scaling of bk with k 
follows a power law. However, for strongly assortative (e.g., 
c = 0.6) networks, the scaling is not as straightforward in 
that the values of b k for large-degree nodes are not necessar¬ 
ily high but instead spread in a relatively wide range. (Simi¬ 
lar phenomena were observed previously in some real-world 
networks [42].) The density of the infected nodes within the 
degree class k is p k = I k /N k , where N k and I k are the num¬ 
bers of nodes and of infected nodes within the degree class, 
respectively. From Fig0b), we see that for dissortative (e.g., 
c = — 0 . 2 ), uncorrelated (c = 0 ), and weakly assortative 
(e.g., c = 0.25) networks, p k increases with k. However, 
for strongly assortative (e.g., c = 0.6) networks, the value of 
p k for large-degree nodes can assume relatively small values. 

Combining Figs. [3|a) and Ob), we deduce that the algo¬ 
rithmic betweenness of a node is positively correlated with 
the probability of its being infected. We define the infection 
probability of a node i as pi = Tj/T, where T is the total 
time lapse and Tj is the number of time steps at which node 
i is infected (here we set T = 1000). In Fig. [4] we plot the 
dependence of p, on bi for different values of c. We see that, 
as a general trend, pi increases with bi and the largest h, cor¬ 
responds to the highest value of pi, regardless of the value of 
c. 

Since the node having the maximum algorithmic between¬ 
ness b max is most likely to be infected, it plays a crucial role 
in determining the traffic-driven epidemic threshold. In partic¬ 
ular, only when the spreading rate is large enough to success¬ 
fully infect the node with b max will the epidemic spreading 
be able to sustain. To confirm this, we calculate the depen¬ 
dence of b m ax on the assortativity coefficient c for networks 
with different values of the degree exponent 7 , as shown in 
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FIG. 5: (Color online) For scale-free networks of different values of 
the degree exponent 7 , the maximum nodal algorithmic betweenness 
bmax versus the assortativity coefficient c. Each data point is result 
of averaging over 100 random network realizations. 



FIG. 6: (Color online) For scale-free networks of degree exponent 
7 = 3 with different values of the assortativity coefficient c, epi¬ 
demic threshold /3 C as a function of the packet-generation rate R. 
The slope of the fitted line is about -1. Each data point results from 
an average over 100 network realizations. 


Fig- 0 We see that, for each value of 7 , bmax peaks at some 
moderate value of c. Comparing with the results in Fig. [3 we 
see that the b max versus c behavior is opposite to that of j3 c 
versus c, implying that increasing b ma x causes /3 C to decrease. 
The critical value of c leading to the maximum b ma x is in 


fact indistinguishable from the value that minimizes (3 C . As a 
by-product, the phenomenon shown in Fig. [3 suggests a rela¬ 
tionship between the maximum algorithmic betweenness and 
the largest eigenvalue of the betweenness matrix, for which a 
mathematical understanding is not yet available. 

Finally, we study the epidemic threshold /3 C as a function 
of the packet-generation rate R for different values of the as¬ 
sortativity coefficient c. The results are shown in Fig. [3 We 
see that /? c scales inversely with R for each value of c, as pre¬ 
dicted by our mean-field theory [Eq. ©I- This is intuitively 
correct, as an increase in the traffic flow generally facilitates 
epidemic outbreak. 


IV. CONCLUSIONS 

To summarize, our investigation of traffic-driven epidemic 
spreading in correlated scale-free networks reveals the exis¬ 
tence of a critical level of network correlation for which the 
outbreak of epidemic is maximally promoted. We use an ex¬ 
tended degree-based mean-field theory, taking into account 
the traffic flow dynamics, to account for the phenomenon. In 
addition, we find that nodes with larger algorithmic between¬ 
ness are more likely to be infected. We note that in tradi¬ 
tional epidemic spreading where infections are transmitted as 
a reactive process from nodes to all neighbors, the epidemic 
threshold is largely determined by the nodes with the largest 
degree |43|]. Our results reveal, however, that for traffic-driven 
epidemic spreading, the threshold mainly depends on the node 
with the maximum algorithmic betweenness. In particular, 
the threshold tends to decrease with the value of the max¬ 
imum algorithmic betweenness. While for uncorrelated or 
disassortative networks, the largest-degree node typically has 
the maximum algorithmic betweenness, this is not the case 
for assortative networks. Our results indicate that in traffic- 
driven spreading dynamics, the structural properties of the 
network can have a significant effect on epidemic spreading 
and outbreak, providing a potential mechanism to control var¬ 
ious spreading dynamics in real-world situations. 
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